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HIGHLIGHTS 


•  Exact  solutions  were  obtained  to  simulate  current  passages  in  SOFC  electrodes. 

•  A  solid  oxide  fuel  cell  (SOFC)  electrode  was  selected  as  a  test  case  for  these  studies. 

•  Assess  geometry  influence  3-D  electrode  microstructure  on  performance. 

•  Performance  was  correlated  to  geometrical  changes  in  the  electrode. 
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Exact  solutions  were  obtained  for  variations  in  the  potential  and  the  current  for  three  axisymmetric 
geometries,  with  positive,  negative  and  zero  curvatures,  which  simulate  current  transport  in  fuel  cell 
electrodes.  These  solutions  can  be  used  to  assess  the  influence  of  geometry  on  performance  for  three 
dimensional  electrode  microstructures.  A  solid  oxide  fuel  cell  (SOFC)  electrode  was  selected  as  a  test  case 
for  these  studies.  From  the  exact  solutions,  simulations  of  current  flow  and  potential  drop  for  one 
dimensional  networks  in  SOFC  electrodes  were  performed.  Numerical  tests  demonstrated  that  surfaces 
with  positive  curvature  have  greater  current  flow  for  the  same  potential  drop  due  to  higher  current 
losses  through  the  lateral  surface  area.  The  study  also  showed  that  zero  curvature  solutions  will  be 
sufficiently  accurate  for  positive  or  negative  curvature  geometries  for  moderate  radius  changes,  but  differ 
significantly  from  positive  or  negative  curvature  solutions  for  more  extreme  radius  changes.  Analytical 
solutions  indicate  fundamental  differences  in  geometry  and  its  influence  on  current  flow.  Based  on  the 
results  of  the  simulations,  an  approximate  solution,  based  on  one  non-dimensional  parameter,  was 
developed  for  estimating  the  effects  of  extreme  changes  in  cross-section  area. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  composite  electrodes  applied  in  electrochemical  energy 
storage  and  conversion  devices  are  heterogeneous  functional 
material  systems  that  must  support  the  effective  movement  of 
charge  generated  by  the  electrochemical  interaction  between 
electrode  phases.  Such  electrodes  may  include  ceramic— metal 
(cermet)  composites  and  mixed  conducting  electrodes  found  in 
solid  oxide  fuel  cells  (SOFCs)  [1],  carbon  and  transition  metal 
oxide  electrodes  of  Li-ion  batteries  [2,3  ,  and  porous  titania  films 
used  in  dye-sensitized  solar  cells  4,5].  These  chemically  active 
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composite  systems  require  the  presence  of  sites  that  support 
fundamental  chemistry  and  sufficient  transport  networks  that 
deliver  reactant  species  to  and  byproducts  from  these  active  sites. 
Considering  these  two  key  requirements,  the  deliberate  design  of 
microstructural  geometry  in  composite  electrodes  can  benefit 
performance  by  increasing  the  efficacy  of  active  sites  and  by 
promoting  the  efficient  movement  of  charge  through  conductive 
phases  to  the  active  surfaces. 

Within  SOFC  electrodes,  the  effective  transport  of  ions  and 
electrons  to  active  sites  is  of  particular  interest.  In  more  common 
cermet  electrodes,  these  sites  reside  where  the  solid  ion  and  elec¬ 
tron  conducting  phases  meet  an  open  pore  phase  that  enables 
transport  of  gaseous  reactants.  In  electrodes  that  support  the  mixed 
conduction  of  ionic  and  electronic  charge  carriers  these  reaction 
sites  may  be  distributed  over  a  broader  portion  of  the  solid-pore 
interface.  It  has  generally  been  found  that  the  extension  of  active 
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Nomenclature 

Pv 

Legendre  function  of  order  v 

R 

charge  transfer  resistance,  Q  m2 

a 

radius  of  cross-section,  m 

R'c 

equivalent  cross-section  resistance,  Q 

Tli,  A2 

arbitrary  constants 

Rs 

equivalent  surface  resistance,  Q 

71c 

area  of  cross-section,  m2 

r0f  re 

entry  and  exit  radii  respectively,  m 

Tleff 

effective  area  of  cross-section,  m2 

X 

transformed  axial  coordinate  /<z  -  z0) 

71s 

surface  area,  m2 

z 

axial  coordinate,  m 

ic 

total  current  through  cross-section,  A 

Zo 

axial  offset,  m 

h  Ki 

modified  Bessel  functions  of  order  one 

a 

hyperbolic  roots,  ±^/l  +  2 /oRk 

k 

curvature  or  slope,  m-1 

P 

equivalent  circuit  coefficient 

L 

segment  length,  m 

<2> 

transformed  potential,  cpja 

£ 

symmetric  half  segment,  length,  m 

potential,  V 

sites  away  from  the  planar  interface  between  electrode  and  elec¬ 
trolyte  layers  of  the  cell  can  reduce  overall  polarization  resistance 
within  the  electrode  and  increase  SOFC  performance.  This  exten¬ 
sion  is  often  achieved  through  the  fabrication  of  a  functional  layer 
near  the  electrode/electrolyte  interface  with  a  finer  microstructure 
than  the  porous  electrode  support.  This  extended  active  surface 
results  in  the  more  effective  transfer  of  ionic  charge  carriers  from 
the  bulk  electrolyte  to  active  sites  within  the  composite  electrode 
[6,7]. 

The  coupled  electrochemistry  and  charge  transport  that  occurs 
in  these  electrodes  have  been  explored  using  a  variety  of  tech¬ 
niques  including  thin-film  electrode  models  [6,8],  extended  surface 
electrode  models  [7,9],  and  percolation  theory  [10-14].  A  com¬ 
monality  in  these  modeling  approaches  is  the  application  of  a 
second  order,  ordinary  differential  equation  comparable  to  the 
equation  describing  heat  transfer  from  a  thermal  fin.  Early  appli¬ 
cation  of  this  equation  can  be  attributed  to  Kenjo  et  al.  [6,8],  who 
used  a  thin-film  model  to  predict  the  polarization  resistance  of 
SOFC  electrodes  with  varying  composition  and  microstructural 
quality.  This  approach  abstracted  the  dispersed  ionic  conducting 
phase  as  a  thin  film  coating  the  electronic  conducting  phase  and 
addressed  both  charge  transport  and  microstructural  geometry 
through  the  use  of  a  single  lumped  parameter.  It  was  found  that 
materials  that  sintered  well,  forming  wide  contact  areas  between 
individual  particles  in  the  composite  electrode,  displayed  superior 
performance  compared  to  materials  with  poor  sinterability.  It  was 
also  demonstrated  that  for  well  sintered  electrodes  polarization 
resistance  could  be  reduced  by  increasing  the  thickness  of  the 
active  electrode.  Flowever,  this  improvement  showed  a  diminishing 
return  that  suggests  a  trade-off  between  transport  losses.  The  thin- 
film  model  of  Kenjo  et  al.  was  adapted  to  an  extended  surface 
formalism  by  Tanner  et  al.  [7],  who  treated  the  ionic  conducting 
phase  as  a  set  of  rectangular  corrugations,  coated  with  electro¬ 
catalyst,  extending  from  the  bulk  electrolyte.  This  approach  added 
microstructural  detail  by  moving  away  from  the  lumped  parameter 
description  of  charge  transport  and  microstructure.  Instead, 
microstructure  was  accounted  for  through  the  sizing  and  separa¬ 
tion  of  the  extended  ion  conducting  surfaces.  Other  descriptions  of 
charge  transport  within  composite  SOFC  electrodes  have  drawn  on 
percolation  theory  to  define  effective  conductivity  and  other  charge 
transfer  parameters  [10-14].  As  in  the  analytical  approaches  pro¬ 
vided  by  Kenjo  et  al.  and  Tanner  et  al.,  these  approaches  apply  a 
governing  equation  comparable  to  that  of  a  constant  cross-section 
thermal  fin. 

The  effects  of  varying  particle  contact  geometry  have  been 
investigated  by  Zhao  and  Virkar  [15  through  the  development  of 
analytical  expressions  for  the  effective  conductivity  of  sintered 
electrodes.  The  incorporation  of  more  detailed  microstructural 
geometry  within  extended  surface  electrode  models  was  also 


explored  by  Nelson  et  al.  [9]  through  the  use  of  a  governing  equa¬ 
tion  that  accounts  for  variable  cross-section  geometry  in  the 
extended  surface.  In  this  analysis,  chains  of  ionic  conducting  par¬ 
ticles  extending  from  the  electrolyte  were  simulated  based  on  an 
iterated  series  of  truncated  cones.  As  illustrated  in  Fig.  1,  these 
particle  chains  can  be  considered  to  comprise  the  active  layer  of  the 
electrode  that  resides  near  the  electrode-electrolyte  interface.  In 
electrodes  with  good  connectivity  such  chains  may  extend  into  the 
porous  electrode  support.  The  variable  cross-section  extended 
surface  approach  was  found  to  be  capable  of  replicating  experi¬ 
mentally  observed  polarization  resistances  in  SOFC  electrodes,  in¬ 
clusive  of  resistive  effects  linked  to  reduced  interparticle  contact 
areas  that  result  when  electrode  particles  are  not  sufficiently  sin¬ 
tered  [9  .  Further  investigation  of  the  impacts  of  multi-dimensional 
transport  and  space  charge  regions  at  contact  points  between 
particles  demonstrated  that  the  analytical  extended  surface  models 
provided  an  accurate  assessment  of  electrode  performance, 
particularly  for  purposes  of  screening  potential  microstructural 
designs  [16].  Thus  it  has  been  found  that  variable  cross-section 
extended  surface  models  present  a  promising  tool  for  the  more 
deterministic  design  of  electrode  microstructures.  However,  the 
development  of  a  more  robust  means  of  assessing  electrode  quality 
can  benefit  from  extension  of  the  electrochemical  fin  modeling 
concept  to  a  more  diverse  set  of  particle  shapes  that  may  be  present 
in  electrode  microstructures.  For  example,  the  fabrication  of  elec¬ 
trodes  based  on  spherical  particles  has  been  demonstrated  using 
samaria-doped  ceria,  and  the  conductivity  of  these  porous  com¬ 
posites  has  shown  clear  sensitivity  to  microstructural  geometry 
[15]. 

In  the  present  work  exact  analytical  solutions  are  presented 
for  a  variety  of  particle  shapes  that  may  be  used  to  simulate  an 
extended  surface  electrode.  These  shapes  include  conical  particles 
similar  to  those  investigated  previously  [9,16],  as  well  as  particles 
with  spherical  shapes  and  particles  with  a  profile  based  on  a 
hyperbolic  cosine  function.  These  latter  two  particle  profiles  are 
illustrated  in  Fig.  lb.  Using  a  set  of  basic  parametric  studies, 
transport  effects  associated  with  particle  curvature,  variable 
cross-sectional  area,  surface  area,  and  length  are  investigated  to 
further  elucidate  the  role  of  microstructural  geometry  in  elec¬ 
trode  performance.  The  effects  of  material  properties  and 
boundary  conditions  are  included  in  the  solutions  and  can  also  be 
addressed  to  better  determine  the  interaction  between  particle 
shape  and  the  physical  phenomena  that  govern  electrode  oper¬ 
ation.  For  example,  the  effects  of  curvature  on  surface  current 
losses  can  be  examined  since  exact  solutions  for  positive 
(spherical  surface),  zero  (conical  surface)  and  negative  (hyper¬ 
bolic  surface)  curvatures  have  been  developed.  Hence,  even  small 
differences  in  current  loss  can  be  used  to  find  trends.  The 
analytical  solutions  have  shown  conclusively  that  positive 
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Fig.  1.  To  formulate  an  extended  surface  electrode  model,  the  active  layer  of  a  composite  electrode  (a.)  is  treated  as  a  network  of  particle  chains  supporting  the  conduction  of  ions 
away  from  the  bulk  electrolyte.  The  particles  comprising  such  chains  may  be  approximated  using  several  shapes  including:  conical  (a.),  spherical  (a.  and  b.),  and  hyperbolic  (b.).  The 
governing  equation  for  charge  transport  in  these  variable  cross-section  structures  may  be  developed  by  taking  a  charge  balance  on  a  differential  control  volume  (c.)  within  the 
structure. 


curvatures  will  have  a  greater  loss  than  negative,  with  zero  cur¬ 
vature  producing  an  intermediate  result. 

2.  Model  formulation 

The  governing  equations,  that  establish  the  analytical  model 
presented  here,  are  developed  based  on  an  extended  surface  of 
charge  conducting  material  subject  to  a  charge  transfer  reaction  at 
its  surface.  As  illustrated  in  Fig.  1,  the  profile  of  this  extended 
surface  may  take  a  number  of  shapes.  However,  to  permit  the 
cross-sectional  averaging  of  the  potential  within  the  structure, 
these  shapes  must  be  axisymmetric.  Two  key  assumptions  are 
applied  with  respect  to  transport  in  the  solid  cross-section  of  the 
extended  surfaces  considered.  First,  one  dimensional  transport  is 
assumed  along  the  length  of  these  extended  surfaces  in  the  z-di- 
rection.  In  previous  investigations  this  approximation  has  been 
found  to  hold  for  most  cross-sectional  geometries  [16].  An 
exception  to  the  sufficiency  of  this  assumption  is  found  for  par¬ 
ticles  with  contact  radii  that  are  significantly  smaller  than  the 
largest  particle  radius,  i.e.  a  high  gradient  in  particle  cross- 
sectional  area.  In  such  cases,  the  one-dimensional  approximation 
overestimates  performance  predictions.  Second,  in  the  analytical 
model  presented  the  flow  of  charge  within  the  solid  phase  is 
treated  as  ohmic,  with  space  charge  effects  neglected.  In  contact 
regions  between  individual  particles,  grain  boundaries  may  lead  to 
the  formation  space  charge  regions,  which  display  a  higher 
resistance  to  transport  than  the  material  in  the  bulk  particle 
[15,17].  As  in  the  case  of  multidimensional  transport,  neglecting 
the  additional  resistance  introduced  by  these  space  charge  regions 
leads  to  an  overestimate  of  electrode  performance  by  the  analyt¬ 
ical  models  [16].  In  general,  the  analytical  approach  presented  in 


this  work  proves  sufficient  as  a  meso-scale  screening  tool  for 
assessing  the  quality  of  SOFC  electrode  microstructures  that  may 
be  produced  using  emerging  fabrication  techniques. 

Considering  the  case  of  one-dimensional  transport  in  the 
absence  of  space  charge  effects,  the  flow  of  current  within  the  solid 
cross-section  may  be  defined  according  to  Ohm’s  law. 

ic  =  -<rAc(z)^  (1) 

Here,  a  is  the  conductivity  of  the  solid  phase,  Ac  is  the  cross- 
sectional  area,  and  cp  is  a  potential  difference  between  the  solid 
phase  and  a  reference  potential.  Details  of  the  definition  of  this 
reference  potential  are  provided  by  Nelson  et  al.  [9]. 

To  permit  development  of  an  analytical  model  a  linear 
description  of  the  surface  charge  transfer  reaction  is  applied.  For  an 
SOFC  electrode  this  approximation  is  considered  to  hold  when  gas 
transport  is  fast  within  the  pores  of  the  electrode.  In  such  a  case,  the 
relative  rate  of  charge  transfer  may  be  approximated  using  a  charge 
transfer  resistance.  Based  on  this  approximation  a  current  related 
to  the  surface  charge  transfer  may  be  defined  in  terms  of  the  po¬ 
tential  difference,  cp,  the  surface  area,  As,  and  the  charge  transfer 
resistance,  R. 


As 

With  the  solid  phase  current,  ic,  and  surface  charge  transfer 
current,  is,  defined  the  governing  equation  for  the  analytical  model 
may  be  developed.  This  is  done  by  defining  a  differential  control 
volume  within  the  extended  surface  structure  and  subsequently 
taking  a  charge  balance  on  this  volume,  as  illustrated  in  Fig.  1  c.  For 
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this  development  the  cross-sectional  and  surface  areas  are  defined 
for  a  generic  axisymmetric  structure  according  to  the  local  radius, 
a(z),  as  shown  in  Eqs.  (3)  and  (4).  The  charge  balance,  in  terms  of 
the  currents  is  given  according  to  Eq.  (5). 


2d  cp  d cp  2 


X^-r-^  +  X-^-X^Cp  =  0. 

dx2  dx 


The  general  solution  of  Eq.  (10)  is 


(10) 


Ac(z)  =  7 ra2(z) 


dAs(z)  =  2na(z)\l  1  +  (/J“j  dz 


*C  +  j/4z  )  -  i'c  - 


2  na(z)yjl  +  (^)2d  z 


<p  =  0 


(3) 

(4) 

(5) 


<P  =  A^Ii  (x)  +A2Kj  (x) 


=  APi{2\/^jP-(kz  +  z0) 


+  AlI9  I  -(kz  +  z0) 


HD 


Substituting  Eq.  (1)  for  the  solid  phase  current  and  rearranging 
terms  results  in  an  ordinary  differential  equation  describing  the 
potential  difference  within  the  structure. 


where  A\  and  A2  are  arbitrary  constants,  evaluated  by  the  boundary 
conditions  at  x  =  0,  and  x  =  L,  and  Ji,  K\  are  modified  Bessel  func¬ 
tions  of  order  one  of  the  first  and  second  kind,  respectively. 


a2(z) 


d  2cp 
dz2 


+  2  a(z) 


da  d cp 
dz  dz 


2a(z) 

oR 


(6) 


In  Eq.  (6)  the  conductivity,  cr,  and  the  charge  transfer  resistance, 
R,  are  constant  functions  of  z.  Eq.  (6)  governs  the  variation  in  the 
potential  for  axisymmetric  geometries  for  a  radius  that  is  an  arbi¬ 
trary  function  of  axial  position.  Note  that  there  is  a  multiple  of  the 
radius  in  each  of  the  three  terms  that  can  be  removed.  Boundary 
conditions,  for  Eq.  (6),  include  a  specified  potential,  rp,  a  specified 
current,  i0  or  a  linear  combination  of  the  potential  and  current. 


3.  Exact  solutions 


3.2.  Spherical  solution 


For  the  case  of  a  spherical  segment  the  radius  is 

a2  +  (z-z0)2=T  (12) 

Substituting  Eq.  (12)  into  Eq.  (6)  the  governing  equation  for  the 
potential  becomes 


d 2  cp 


d  cp 


2  ( Z  Z°)  )  Ar,  O’ 


dz2 


dz  aRk' 


(13) 


Exact  solutions  remove  the  uncertainty  from  numerical  ap¬ 
proximations,  allowing  interpretations  to  be  confidently  made 
from  small  differences  in  the  results.  To  determine  the  influence  of 
various  parameters  geometries  have  been  chosen  for  which  exact 
solutions  can  be  found.  However,  the  square  root  in  Eq.  (6)  makes  it 
difficult  to  find  exact  solutions  for  the  potential.  Exact  solutions  for 
three  geometries,  including  spherical  segments,  conical  segments 
and  a  hyperbolic  geometry,  have  been  developed  as  an  aid  in 
determining  the  influence  of  curvature.  Three  parameters  will  be 
used  to  describe  the  geometry:  1 )  the  length  of  the  segment,  L,  2)  a 
curvature,  or  slope,  parameter,  /<,  and  3)  an  offset,  Zq. 

3.1.  Conical  solution 

Conical  segments  have  a  linear  variation  in  the  radius,  making 
the  square  root  a  constant  in  Eq.  (6)  and,  thus,  eliminating  the 
complication  from  the  square  root.  We  use  for  the  radius 


a  =  kz  +  z0.  (7) 

Transforming  the  independent  variable  from  z  to  a  Eq.  (6) 
becomes 


2cp  ndcp  2\/\  +  /<2 

da2  da  aRk2  ^ 


(8) 


Eq.  (8)  can  be  converted  to  a  standard  Bessel  function  differ¬ 
ential  equation  by  setting 


x  =  2 


2(1  +/<2)) 


aRk 


a. 


(9) 


Eq.  (13)  can  be  simplified  by  the  substitution 


x  =  k(z  -  z0), 


(14) 


which  yields  the  form 
42 


0 

Eq.  (15)  is  the  Legendre  equation  and  has  the  solution 


(15) 


cp  =  A,Pv{r](x)  +A2PV(2)(X) 

=  A^pv^(k(z  -Zq))  A2Pv(2)(k(Z  -Zq))  (16) 


where  A\  and  A2  are  arbitrary  constants  and 

p(l ,  2)  =  -1/2  ±  yj (1/2)2  +  (2 /aRk).  A  series  solution  about  x  =  0 
was  developed  for  this  equation.  One  series  begins  at  x°  and  con¬ 
tains  only  even  powers  of  x.  The  second  begins  at  x1  and  contains 
only  odd  powers  of  x.  The  radius  of  convergence  of  the  series  is  one 
and  hence  does  not  extend  past  the  singular  points  x  =  ±1.  That  is, 
the  series  does  not  converge  at  or  beyond  the  poles  of  the  sphere. 
This  convergence  behavior  is  considered  sufficient  because  a  par¬ 
ticle  that  is  a  complete  (i.e.  non-truncated)  sphere  would  represent 
an  unsintered  geometry,  and  therefore  be  a  non-functional  elec¬ 
trode  microstructure. 


3.3.  Hyperbolic  solution 

A  solution  for  a  hyperbolic  geometry  can  be  attained  by  noting 
that  if  da/dz  =  sinh(/<(z  -  a))  (i.e.,  take  a(z)  =  1 //<cosh(/<(z  -  z0))) 

then  yj  1  +  (da/dz)2  =  cosh(/<(z  -  z0)).  For  cp  proportional  to 
sinh(/<z  -  z0))  each  term  in  Eq.  (6)  will  contain  the  factor 
cosh2(/c(z  -  Zo))sinh(/<(z  -  Zq)). 


Eq.  (8)  becomes 
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Using 


a(z)  =  ^cosh  (fc(z-zo)) 


(18) 


and  substituting  Eq.  (18)  into  Eq.  (6)  the  governing  equation  for  the 
potential  is 


('+*2)w+3*,3rm’’-0 


dcp 


(19) 


where  x  =  sinh(/<z  -  zo)).  Note  the  similarity  between  Eq.  (19)  and 
Eq.  (15).  The  change  in  sign  of  x2  in  the  first  term  illustrates  the 
change  in  curvature  of  the  surface.  As  in  the  case  of  the  spherical 
segment,  a  series  solution  about  x  =  0  can  be  developed.  One  series 
will  begin  at  x°  and  contain  only  even  powers  of  x.  The  second  will 
begin  at  x1  and  contain  only  odd  powers  of  x.  Unfortunately,  the 
radius  of  convergence  of  the  series  is  one  and  hence  does  not 
extend  past  x  =  ±1.  The  poles  in  this  case  are  at  x  =  ±i,  but  the  series 
diverges  when  |x|  =  |sinh(k(z-z0))|  >  1.  Hence  a  more  general 
solution  needs  to  be  developed  to  move  beyond  x  =  ±1 

A  new  solution  can  be  developed  by  changing  the  dependent 
variable  in  Eq.  (6).  In  a  manner  analogous  to  the  solution  of  the 
transient  diffusion  equation  in  spherical  coordinates,  let 


0 


cp  = 


a(z) 


Eq.  (6)  becomes 


(20) 


d2<f>  1 


da  2  < 


S' 


0  =  0. 


dz2  a 

Substituting  for  the  radius  from  Eq.  (18)  results  in 


d2<Z> 

dz2 


- k 2 


1 


2 

aRk 


0  =  0. 


This  equation  can  be  readily  solved  to  yield 

0  =  A1ea*+A2e-a*, 


(21) 


(22) 


(23) 


where  a  =  sj\  +  2 /oRk,  x  =  k(z  -  zo),  and  A\,  A2  are  arbitrary 
constants. 


3.4.  Composite  solution 

The  individual  general  solutions,  Eqs.  (11),  (16)  and  (23),  can 
readily  be  combined  to  build  a  composite  structure.  Each  segment 
of  the  structure  is  described  by  the  nodes  at  each  end,  the  geo¬ 
metric  parameters  ( particle  length,  particle  curvature  and  offset ),  the 
conductivity,  and  the  charge  transfer  resistance.  The  potential  at 
the  nodes  are  the  unknowns,  and  the  governing  equations  are 
generated  by  requiring  continuity  of  the  current  (i.e.,  the  sum  of  the 
currents  entering  a  node  minus  the  sum  of  the  currents  leaving  a 
node  is  zero.) 


4.  Parametric  simulations 

Computer  codes  have  been  written  based  on  these  solutions  to 
simulate  arbitrary  structures  using  the  exact  solutions  for  a  conical, 
spherical  and  hyperbolic  geometry.  Two  versions  of  the  code  have 
been  developed:  1 )  a  MATLAB  version  and  2)  an  Excel  version  offset 
with  Visual  Basic.  Each  version  has  its  own  advantages  and  disad¬ 
vantages.  The  Excel  version  is  easily  ported  from  machine  to 


machine.  The  MATLAB  version  can  readily  accommodate  a  model 
with  many  segments. 

Nine  cases  with  realistic  parameters  were  considered  as  shown 
in  Fig.  2.  There  are  three  cases  for  each  of  the  signs  in  the  curvature. 
In  each  case  the  radius  varies  by  a  factor  of  two  from  the  minimum 
to  the  maximum.  The  first  case  consists  of  one  segment  with  the 
larger  radius  at  twice  the  potential  of  the  smaller  radius.  The  cur¬ 
vature  geometry  also  has  two,  two  segment  versions:  one  with  the 
minimum  radius  at  the  center  and  one  with  the  maximum  radius  at 
the  center. 

Fig.  2  shows  that  the  input  current  is  approximately  the  same  for 
each  of  the  two  segment  cases,  but  there  is  a  significant  difference 
for  the  single  segment  cases.  This  result  appears  to  be  due  to  the  fact 
that  the  average  area  of  the  cross-section  for  the  spherical  case  is 
larger  while  the  hyperbolic  case  has  a  smaller  average  area.  The 
conical  segment  is  intermediate.  The  current  loss  out  the  surface 
though  is  largest  for  the  spherical  segment,  smallest  for  the  hyper¬ 
bolic  segment,  and  again  is  intermediate  for  the  conical  segments. 
The  two  segment  cases,  when  compared  for  each  geometry,  show  it 
makes  little  difference  where  the  minimum  radius  is  located. 

Fig.  3  shows  an  investigation  of  potential  as  a  function  of  axial 
location  in  the  three  geometries:  spherical,  conical,  and  hyperbolic. 
It  appears  that  the  conical  segment  provides  a  good  approximation 
over  moderate  changes  in  radius,  and  should  provide  accurate 
current  predictions  for  the  spherical  and  hyperbolic  cases  if  enough 
conical  segments  are  used  across  a  single  spherical  or  hyperbolic 
segment  to  capture  the  non-zero  curvature. 

To  evaluate  how  neck  size  affects  the  influence  of  geometry  on 
the  current  flow,  a  radius  ratio  parameter  is  defined  as  the  ratio  of 
the  end  radius  of  a  particle  to  the  neck  radius  that  serves  as  the 
contact  area  between  two  artificially  sintered  particles.  Varying  this 
radius  ratio  about  that  which  is  shown  in  Fig.  2  for  each  of  the  two- 
segment  cases  allows  the  model  to  assess  the  influence  that 
constrictive  effects  have  on  current  flow  for  these  different 
geometries. 

Fig.  4  is  a  plot  of  resistance  per  unit  length  that  results  from 
setting  a  constant  current  flux  boundary  condition  at  one  end  of  the 
two-segment  particles  and  fixing  a  zero  potential  at  the  opposite 
end.  This  value  is  plotted  for  varying  radius  ratios  for  when  the 
minimum  radius  and  maximum  radius  are  at  the  center,  respec¬ 
tively.  In  this  figure  it  appears  that  increased  deviation  from  the 
cylindrical  case  of  rend/rneck  =  1  yields  a  higher  equilibrium  po¬ 
tential,  regardless  of  the  curvature  geometry.  This  is  due  to 
constrictive  effects  of  the  reduced  neck  size  between  sintered 
particles  when  the  minimum  radius  is  at  the  center,  and  reduced 
cross  sectional  area  for  current  to  flow  at  the  ends  of  the  sintered 
particles  when  the  maximum  radius  is  at  the  center. 

The  results  in  Fig.  4  also  indicate  that  as  the  radius  ratio  deviates 
further  from  the  value  of  two,  presented  in  Fig.  2,  the  aforemen¬ 
tioned  differences  between  curvature  geometries  are  more  pro¬ 
nounced.  This  becomes  most  apparent  for  extreme  ratios  of  0.1  and 
10,  where  greater  differences  in  cross  sectional  area  variation  for 
each  geometry  type  more  significantly  affect  the  resulting  current 
flow.  In  these  cases,  hyperbolic  particles  exhibit  the  highest  resis¬ 
tance  per  unit  length,  while  spherical  particles  are  lowest.  Conical 
particles  showcase  intermediate  results.  This  implies  that  along  the 
primary  transport  direction,  spherical  particles  provide  the  least 
resistance  to  the  flow  of  ions  due  to  a  greater  cross  sectional  area 
along  their  length  than  corresponding  conical  and  hyperbolic 
particles. 

5.  Approximate  solution 

The  results  presented  can  be  used  to  provide  insight  into  how 
microstructure  affects  the  performance.  Also,  it  would  be  very 
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Fig.  2.  Current  loss  through  the  surface  of  various  composite  structures  illustrates  that  a  positive  curvature  (spherical  segment)  will  have  a  greater  loss  than  a  zero  (conical  segment) 
or  negative  (hyperbolic  segment)  curvature.  In  all  nine  cases  the  larger  radius  is  twice  the  smaller.  All  segments  have  the  same  length.  The  current  enters  from  the  left  surface  (fn), 
and  exits  to  the  right  surface  (/out).  Current  is  lost  through  the  side  surfaces  as  indicated  by  I\oss.  All  current  values  have  units  of  A  m-2.  One  or  two  particles  were  connected  as 
shown.  The  area  ratio  (maximum  to  minimum)  was  four  in  all  cases.  The  conductivity  (10  S  m-1)  and  surface  loss  coefficient  (10-4  fi)  was  the  same  for  all  examples. 


Position,  x  -  m 


Fig.  3.  The  effect  of  curvature  on  potential  in  a  single  segment  illustrates  that  more  rapid  changes  in  potential  indicate  more  current  loss  through  the  surface.  The  radius  at  current 
entry  (x  =  0)  is  twice  as  large  as  the  radius  of  at  current  exit. 
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Radius  Ratio,  (a  / a  ) 

’  v  p  rr 


Fig.  4.  Resistance  per  unit  length  shown  as  a  function  of  radius  ratio  for  sintered  particles  of  each  geometry  illustrates  that  as  the  radius  ratio  deviates  from  the  case  of  a  cylinder 
(rend/fneck  =  1 ),  current  flow  is  increasingly  impacted  by  the  curvature  geometry. 


useful  for  design  estimates  to  develop  an  accurate  approximate 
analysis.  We  can  develop  an  estimate  by  considering  the  current 
flowing  in  the  circuit  in  Fig.  5.  In  the  figure  current  leaving  the 
surface  is  assumed  to  go  to  the  ground  potential  through  the 
resistor  at  the  center  of  the  segment.  The  conduction  path  with 
currents  in,  7in,  and  out,  Jout,  has  a  total  resistance  of  R'c  =  L/(aAeff) 
where  a  is  the  conductivity,  I  is  the  total  segment  length,  and  Ae ff  is 
the  effective  conduction  cr-sectional  area.  The  resistor  to  ground 
has  a  total  resistance  of  RSIAS,  where  As  is  the  total  surface  area  but 
is  weighted  by  the  potential.  If  the  current,  over  an  average,  is 
assumed  to  see  half  the  potential  then  the  effective  surface 
resistance  is  R's  =  2 Rs/As.  The  circuit  can  now  be  readily  solved  to 
yield 

h  =  l/2[VW)](As/>W)  (24) 

'in  l  +  l/2[I/(ffRs)](As/Aeff) 


where  Is  is  the  current  lost  through  the  surface.  In  Eq.  (24)  we  have 
taken  the  surface  area  to  be  the  actual  surface  area  of  the  segment, 
that  is 


For  the  equivalent  cross-sectional  area  consider  the  current,  /, 
through  a  cross-section  where 


I  =  ff4:(x)g  (26) 

and  >4c(x)  is  the  cross-section  area  at  location  x  at  potential  0. 

The  total  change  in  the  potential  across  the  length  L,  for  no  loss 
through  the  surface  ( Rs  =  oo ),  is 


hn  IOul 


Fig.  5.  Equivalent  resistor  circuit  for  a  single  segment  consists  of  three  resistors.  The 
surface  losses  are  represented  by  the  resistor  Rs  and  the  axial  resistance,  Rc,  is  divided 
into  two  equal  parts  before  and  after  resistor  Rs.  The  currents  in  and  out  of  the  segment 
are  I\n  and  I0ut  respectively,  and  the  current  lost  through  the  surface  is  /s. 


(27) 


which  can  be  used  to  define  the  equivalent  cross-section  area,  Ae ff, 
as: 


1  _  1  f  dx 

Ae ff  L  J  71c  (x) 

o 


(28) 


For  the  three  geometries  for  which  we  have  exact  solutions  the 
effective  cross-section  area  and  the  total  surface  areas  can  be  found 
analytically.  For  the  conical  geometry  symmetric  about  the  origin, 
x  =  0,  the  radius  varies  according  to: 
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r(*)  =  re+(^)x, 


(29a) 


A  -  — 
5  k2 


/<£  +  -sinh(2/<£) 


(31c) 


the  effective  cross-section  over  -C  <  x  <£  is  given  by: 

Aeff  =  7T  r0re,  (29b) 

and  the  surface  area  is 


As  =2^1  +  Vo  +  re)fi.  (29c) 

For  the  spherical  geometry  symmetric  about  the  origin  the 
radius  varies  according 

r(x)  =  VR 2  -x2,  (30a) 

the  effective  cross-section  over  -£  <x  <£  is  given  by: 


_  27T R8. 

B) 

and  the  surface  area  is 


(30b) 


As  =  4  nRi. 


(30c) 


Last  for  the  hyperbolic  geometry  symmetric  about  the  origin  the 
radius  varies  according  to: 


r(x)  =  icosh  (lex), 


(31a) 


Note  that  for  all  three  cases  the  total  length  L  =  2£  Eqs.  (29)— (31 ) 
can  be  used  to  arrive  at  values  for  combinations  of  geometries  by 
applying  the  definition  for  the  effective  cross-section  area. 

Eq.  (24)  can  now  be  generalized  by  using  a  parameter,  (3,  to 
absorb  some  of  the  inaccuracies  due  to  the  assumption  that  the 
surface  resistance  is  connected  at  the  center  of  a  symmetric 
segment.  The  approximation  can  be  written  as: 

k  =  (32) 

hn  1  +  j  2[L/ (cjRs)}(As  / Aeffj 

A  least  square  error  fit  was  used  to  find  from  simulations  on 
the  symmetric  geometries  presented.  The  best  fit  was  given  by  a 
value  of  1.098  for  /?.  The  best  fit  for  the  symmetric  geometries 
considered  is  shown  in  Fig.  6.  The  results  in  Fig.  6  appear  accurate 
and  are  qualitatively  the  same  as  those  presented  in  Fig.  4  of 
reference  [18].  Flence  the  approximation  in  Eq.  (32)  can  be  used  to 
accurately  predict  the  fin  effectiveness.  A  single  parameter,  [1/ 
(<7fls)](As/Aeff),  appears  to  be  sufficient  to  represent  the  response  of 
symmetric  geometries.  The  parameter  is  the  product  of  the  axial 
resistance  relative  to  the  surface  resistance  and  the  surface  area 
relative  to  the  axial  area.  Values  for  the  parameter  only  need  to  be 
accurate  to  within  a  factor  of  1.7.  A  value  of  1.098  for  (3  represents  a 
resistance  of  0.551^  in  Fig.  5.  Note  that  if  (3  is  exactly  one  then  the 
approximation  in  Eq.  (24)  would  result,  and  the  axial  resistance 
would  connect  to  the  surface  resistor  at  its  center  as  in  Fig.  5. 


and  the  effective  cross-section  over  -£  <x  <£  is  given  by: 


^eff  — 


7TC 

k  tanh(fcfi)  ’ 


(31b) 


and  the  surface  area  is 


6.  Conclusions 

Exact  solutions  have  been  obtained  for  the  variation  in  the  po¬ 
tential  and  current  for  three  axisymmetric  geometries:  a  spherical 
geometry  (positive  curvature),  a  conical  geometry  (zero  curvature) 


Fig.  6.  A  comparison  of  the  approximate  solution  using  /?  =  1.098  with  the  results  for  various  symmetric  geometries  accurately  predicts  the  fin  effectiveness  for  the  symmetric 
geometries  investigated.  Note  that  most  of  the  performance  change  occurs  over  a  single  decade  in  the  nondimensional  parameter  0.25  <  Bi(As/Ac)  <  2.5  which  agrees  with  the 
results  presented  in  Ref.  [18]. 
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and  a  hyperbolic  geometry  (negative  curvature).  The  solutions  can 
be  used  to  accurately  assess  the  influence  of  geometry.  Geometries 
with  positive  curvature  have  more  current  flow  for  the  same  po¬ 
tential  drop  and  the  same  change  in  radius,  but  lose  more  current 
through  the  lateral  surface  area.  Conical  segments  will  be  suffi¬ 
ciently  accurate  for  moderate  changes  in  radius  if  several  segments 
are  used  to  represent  a  single  spherical  or  hyperbolic  segment.  The 
mathematics  in  the  derivation  of  the  solutions  illustrate  the 
fundamental  difference  in  the  geometries  and  the  influence  of  ge¬ 
ometry  on  potential  variations  and  current  flow. 

An  approximate  solution  was  developed  that  provides  a  quick 
and  accurate  assessment  of  changes  in  geometry  and  material 
properties.  A  single  dimensionless  parameter  combining  physical 
parameters  and  geometric  parameters  appears  to  be  sufficient  for 
representing  the  performance  of  symmetric  particles.  This  can 
provide  insight  into  the  effects  of  microstructure  on  performance 
and  can  aid  in  the  design  of  new  microstructures  without  the  need 
for  detailed  analysis. 
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